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Abstract 

The propagation of an acoustic wave through two-phase porous media with spatial variation 
in porosity is studied. The evolutionary wave equation is derived, and the propagation of 
an acoustic wave is numerically analyzed in application to marine sediments with different 
physical parameters. 



Introduction 


Investigation of acoustic-wave propagation in two-phase porous media, marine sediments in 
particular, finds an increasing interest in studying physical properties of such media. There 
has been a significant amount of publications on the propagation of sound in the sea floor. 
Acoustic-wave propagation in sediments is controlled by intrinsic properties of a sediment, 
characterized by a number of physical parameters. One of the parameters that significantly 
influences sound propagation in a sediment is the porosity. The porosity indicates relative 
amounts of solid and liquid fractions in a sediment and hence it determines the frame bulk 
and shear moduli and through this the acoustic-wave speed. The variations in the medium 
properties can arise due to random packing of inhomogeneous sediment grains. In Refs. 
[DIMS] dispersion and acoustic-wave scattering from randomly varying heterogeneities in 
the poroelastic medium properties such as the porosity or the frame bulk modulus has been 
studied experimentally and theoretically. In Ref. [4J perturbation theory is used to derive 
a poroelastic wave equation which describes the first-order scattering by the heterogeneities 
of a medium. The scattering of energy from heterogeneities accounts for additional losses 
of a sound wave propagating through a poroelastic medium. It was shown that random 
variations in the parameters have more significant effects on the sound propagation through 
a consolidated medium than through a sand sediment. 

In the present paper, we study the influence of porosity variations in space on the prop¬ 
agation of a plain acoustic wave in marine sediments. We shall describe the porosity of a 
sediment and similarly some other physical parameters as fluctuations about their average 
values. If the fluctuations are small the sound field is assumed to change slowly at the 
wave-length scale, and we shall use the method of slowly varying form of the wave (see 
Ref. [5]) to develop a poroelastic evolutionary wave equation in a heterogeneous medium. 
This method is similar to the method of a slowly varying amplitude widely used in studying 
nonlinear wave interactions in media with dispersion, for instance, in studying light-wave 
interactions in nonlinear optics. To analyze different cases of spatial variation in porosity 
computer simulation of the obtained evolutionary equation is performed. 


1. Derivation of approximate wave equations 
for finite-amplitude acoustic waves 
in marine sediments 

To derive the approximate evolution wave equation we start from the continuity equations 
for the densities and momenta of the liquid and solid phases of a sediment composed of 
a rigid frame and pores filled with water, see Ref. [6], These equations are equivalent in 
the main features to the equations developed by Biot mm, but they are presented in a 
somewhat different form, we shall not list them here. On the basis of these equations, in Ref. 
[TO j the equations for the densities of the liquid and solid phases, pf and p s , were derived (in 
this paper we don’t take diffraction and nonlinearity into account): 
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In these equations the following variables are used, 


x = ex 
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and the moving coordinate 


r = t — x/c. 


(4) 


(In Eqs. (JI]), ([2]) and everywhere below the primes for x are omitted.) In the relation (J3J) 
the small parameter e is introduced as 


e ~ v x /c ~ u x /c ~ 5pf/p f ~ 5p s /p s , 


(5) 


here c is the speed of sound in the sediment and v , u are the hydrodynamic velocities of the 
liquid and solid phases, 5pj, 5p s are the deviations from equilibrium values of the densities 
of the liquid and solid phases. In Eqs. ©, a the left-hand sides are of the order of ~ e, 
the right-hand side terms are of the order of ~ e 2 . The introduction of the new variables 
©, dH) actually signifies the application of the method of slowly varying wave profile. 

In the equations dU), (J2J) we wrote p n , p s instead of 5p n ,5p s , m is the porosity; 
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where kf, k s and k are the bulk moduli of the fluid, mineral grains constituting the frame, 
and of the frame itself; /i is the shear modulus of the frame; v = 1 — m — k/k s . 

Let us now describe the porosity of a sediment as fluctuations about an average value, 
m = mo + Am. Alongside the porosity the bulk and shear moduli of the frame and the 
speed of sound should similarly fluctuate about their average values, k = k 0 — A k, p = 
p 0 — A/i, c = c 0 — Ac. The deviations of the moduli and speed of sound from their 
equilibrium values have an opposite sign as to the fluctuations of the porosity since with the 
increase of the porosity the sediment frame becomes softer. 

Let us eliminate one of the variables, Spj or 5p s , from the left-hand sides of Eqs. dU), d2])), 
(let it be, e.g., Sp s ) by subtracting one equation from the other one. In the right-hand sides 
the quantity 5p s is expressed through Spf with the formula which is valid to an accuracy 
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Note, that Eqs. (GQ), (J2]) allow two independent longitudinal modes, the so called fast and 
slow waves. As it is shown in Ref. HP, the slow wave (unlike the fast one) is a strongly 
attenuated diffusion mode, and it does not contribute significantly to the sound held. In this 
approximation we arrive at the equation for an acoustic wave in a sediment with parameters, 
that vary with distance: 
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To obtain these equations we took into account that in sand sediments the bulk modulus 
k s of quartz grains is much greater than that of the pore water, and in this case G can 
be evaluated as G ~ m/kf, provided m is not close to zero. In Eq. (J7|) we have used the 
quantity 8(x) which absolute value would be around 0.1, 
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In Eq. (J7J) the term a>iD T pf that accounts for dissipation is introduced. D r is the dissipation 
linear operator in the variable r which is characterized by the property 
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where a is real and positive and it has the meaning of an amplitude attenuation coefficient 
if the coefficient a\ is taken equal to the coefficient at dpf/dx. The relation (J2J) defines 
the action of this operator on any function of the variable r which can be represented by 
a Fourier series or integral. An algebraic expression for a (a;) is a combination of physical 
parameters (complex bulk and shear frame moduli included) of a sediment, and it includes 
the frequency correction function introduced by Biot pjj. 

In Introduction it was noted that acoustic wave scattering from randomly varying hetero¬ 
geneities in the poroelastic medium properties manifests itself in the increase of the sound 
held attenuation. In this paper, we shall not consider sound scattering from heterogeneities, 
including it implicitly in the dissipation term. 


2. Computer simulation of acoustic wave 
propagation in heterogeneous absorbing 
marine sediments 

To study numerically the propagation of acoustic waves in absorbing marine sediments with 
parameters varying with distance we consider Eq. (J7]) presented in a concise form: 

ai ~Bx 4" a i D T pf = 0. (10) 

As it is seen from Eq. (J7D the coefficients ai, ci 2 are the algebraic combinations of physical 
parameters of a sediment. It is convenient to divide Eq. m by ai, 
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Let the density boundary value be 


0. 


pU=o = Ap 0 , A = 10- 3 -10- 5 . 


Introducing a new variable 6 = 10 4 r we have p(27rl0 4 r) = p(2rc9). 

Solving Eq. (HID we are to find the function p(2nQ) periodic in the variable 6 with the 
period 1. The boundary condition is taken to be a harmonic function 


Po = — sin(27T$). 


In Eq. m it is convenient to normalize the functions and the variables except the 
variable x measured in centimeters. We obtain the equation 

| + OT(,)| + B .p = 0. (12) 

p\x=o = -sin(27r6>), (13) 

where 


C = --10%, £ = ±1, 6(x) = 8oS'(x), 

a\ c 

5'e(x) E [0,1], ho ^ [0.1,0.2], 


D e = 10 4 D t 

Eq. (fT2l) describes in fact two processes, the change of the wave phase and the wave 
dissipation. To solve it the so called splitting method P2) is applied. 

Consider a simple example, 


du . 

— = Au + nu, u \ x= o = Uq 

dx 

and calculate u(h), were h is the step in x. 

We can divide the problem into two parts 
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For smooth solutions the equality w(h) = u(h ) + 0[h 2 ) holds true. So, we can obtain the 
solution to Eq. (fT4|) solving two more simple problems. 

The solution of the equation 
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satisfying periodic boundary condition can be received with the difference schemes of the 
’’angle” type. The stencil of this difference scheme is defined by the characteristic equation 
d6 

— = C8'(x). Since S'(x) is positive or equal to zero the direction of the characteristics 
dx 

depends on the sign of the coefficient C. If C > 0 the ’’right angle” difference scheme is 
stable, if C < 0 the ’’left angle” difference scheme is also stable. The condition binding the 

steps in x and 6 is of the form h < —— (h and z are the steps for x and 9 axes respectively). 

\C\ 

If p| x=0 is represented as the Fourier series p| x=0 = Y I/ m(0)e 27r * mei , Eq. ([13]) can be solved 
in an explicit form: 


p(x,6) = ^2 u m(x)e 
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v m {x) = u m ( 0 ) e -“ c /o^'m. 


Since m = ±1 we have p(x,9) = — sin(27r(0 — p)), p — C I 8'(£)d£. This means that the 

Jo 

solution of (1131) gives a shift of the phase equal to 2np. The phase shift moves to the right 
if C is positive and to the left if C is negative. 

Let us consider the equation 

(h + D T p = 0. 


dx 


D t is the linear dissipation operator : 
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10 4 c/ = a is the attenuation coefficient. 

This relation defines the action of this operator on a function of the variable 9 represented 
by the Fourier series p — Y v m e 2mme ■ As a result we obtain the equation 
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p(x,6) = -sm(2iT9)e- a ' lo4x . 
is the propagation distance. 


from which one gets 
Since m = ±1, we have 

The parameter L^ = 

Consider some examples describing the transformation of harmonic acoustic waves prop¬ 
agating along x. It will be seen that the change of the porosity with distance leads to the 
phase shift of the initial acoustic wave. 

In Figures 1 and 2 the graphs of the functions p\ x =o and 
8\x) = 0.667 + 0.333 sin(10” 1 7ra:) are presented. 

We take C = 0.652 • 10^ 2 and 10 4 c/ = a ~ 0.45 • 10~ 2 cm _1 (see experimental data in 
Refs. [131 HT] .IB]). 

The number of nodes of the variable 9 in the interval [0,1] is 64. This number of nodes 
is sufficient for approximating a harmonic function. The step of the spatial variable x is 0.5, 
that corresponds to the magnitude of |C| and the stability condition. 
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The graphs of p(xi, 6) at x t = 25, 50, 75,100 are presented in Figs 3-6. A positive value 
of C gives the wave phase shift to the right. If C is negative (with the same module), the 
wave phase shifts to the left. 

In the above examples the porosity has an oscillating character (Fig. 2) and satisfies the 
inequality 0.331 < S r (x) <1. If C is positive the curve shifts to the right more quickly for 
larger S'(x) values and slowly for smaller S'(x ) values. Negative C values lead to shifting the 
initial curve to the left in a similar way. The porosity as an oscillating function has been 
chosen as an example. 

Now we shall consider more realistic cases of an irregular spatial variation in a sediment 
porosity. Take for example arbitrary continuous 5'{x) functions of the form presented in 
Figs. 7-9. 

Fig. 7: let C > 0. In Figs. 7.1-7.6 the function p(xi,9) is presented for six Xi values: 
Xi= 30, 40, 50, 60, 70 and 80. If 0 < Xi < 40 the initial curve p(xi,9 ) does not shift, that is 
its phase does not change; if 40 < x^ < 60 it moves to the right with an increasing speed; if 
60 < Xi < 100 it moves with a constant high speed. 

Fig. 8. The grafs for the initial function p(xi,9) are not listed here, since they are in 
a sence similar to those in Figs. 7.1-7.6. At 0 < Xi < 40 the phase changes so that the 
function p(xi,9) shifts to the right with a constant speed. At X\ > 40 the speed decreases 
and at Xi > 60 it goes to zero. 

Fig. 9: This case qualitatively repeats the cases of Figs. 7-8 where the curves go respec¬ 
tively up or down. The speed of the initial curve shift slows down at 0 < Xi < 50 and at 
Xi > 50 the shift moves with an accelerating speed. 

In conclusion consider the case of a random pore size distribution (Fig. 10) that can 
arise, for instance, due to random packing of the sediment nonuniform grains. Such porosity 
distribution can be presented as a random digital array, and we shall interpolate it with a 
continuous function, see Fig. 11. This continuous function in its turn can be approximately 
considered as a series of curves of the types presented in Figs. 7-8. That is, the consideration 
given above for Figs. 7-8 can be applied to each section of this function. 

3. Conclusion 

Evolutionary wave equation to describe acoustic wave propagation in a two-phase porous 
media with spatial variations in porosty is derived. Computer simulation of the obtained 
equation is performed to analyze diverse cases of the porosity variations that lead to phase 
shifts in the initial harmonic acoustic wave. 
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FIGURES 



Fig.5. The function p(xi,6 ) at Xi = 75 Fig.6. The function p(xi,9 ) at x t 

















